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We search for the neutrinoless, lepton- flavor- violating tau decays t~ — > £~V°, where I is an 
electron or muon and V° is a vector meson reconstructed as (f> — » K + K~ , p — » ir + ir~,K* — > 
K + tt~ , K — > K~w + . The analysis has been performed using 451 fb _1 of data collected at an e + e~ 
center-of-mass energy near 10.58 GeV with the BABAR detector at the PEP-II storage rings. The 
number of events found in the data is compatible with the background expectation, and upper limits 
on the branching fractions are set in the range (2.6 - 19) x 10" 8 at the 90% confidence level. 

PACS numbers: 13.35.Dx, 14.60.Fg, 11.30.Hv 



Lepton-flavor violation (LFV) involving tau leptons 
has never been observed, and recent experimental re- 
sults have placed stringent limits on the branching frac- 



tions for 2- and 3-body neutrinoless tau decays [3, 0, 0] . 
Many descriptions of physics beyond the Standard Model 
(SM) predict such decays [1, [B| ; and certain models @, Q 
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specifically predict semileptonic tau decays such as r — > 
£ cf>/p/K*/K* (t~ -> ^ _ V°), with rates as high as the 
current experimental limits [3J] . An observation of these 
decays would be a clear signature of physics beyond the 
SM, while improved limits will further constrain models 
of new physics. 

This paper presents a search for LFV in a set of eight 
neutrinoless decay modes r~ — > £~V°{$1, where £ is an 
electron or muon and V° is a neutral vector meson decay- 
ing to two charged hadrons (V° — > h + h~) via one of the 
following four decay modes: <p — > K + K~, p — ► 7r + 7r~, 
K* — > X + 7r~, if — > 7r + if - . This analysis is based 
on data recorded by the BABAR detector at the PEP-II 
asymmetric-energy e + e~ storage rings operated at the 
SLAC National Accelerator Laboratory. The BABAR de- 
tector is described in detail in Ref. ] 
consists of 410 fb -1 recorded at an e + e~ 
(cm.) energy •/& = 10.58 GeV, and 40.. 



The data sample 
center-of-mass 
3fb _1 recorded 
at y/s = 10.54 GeV. With a calculated cross section for 
tau pairs of a TT = 0.919 ± 0.003 nb [l, OH at the stated 
luminosity-weighted y/s, this data set corresponds to the 
production of about 830 million tau decays. 

We use a Monte Carlo (MC) simulation of lepton- 
flavor-violating tau decays to optimize the search. Tau- 
pair events including higher-order radiative corrections 
are generated using KK2f One tau decays via two- 

body phase space to a lepton and a vector meson, with 
the meson decaying according to the measured branch- 
ing fractions [l2| . The other tau decays via SM processes 
simulated with TAU0LA 13]. Final state radiative effects 
are simulated for all decays using PHOTOS 14(. The de- 
tector response is modeled with GEANT4 [HI , and the sim- 
ulated events are then reconstructed in the same manner 
as data. SM background processes are modeled with a 
similar software framework. 

We search for the signal decay t~ — > £~V° — > £~h + h~ 
by reconstructing e + e~ — > t + t~ candidates in which 
three charged particles, each identified as the appropriate 
lepton or hadron, have an invariant mass and energy close 
to that of the parent tau lepton. Candidate signal events 
are first required to have a "3-1 topology," where one 
tau decay yields three charged particles, while the sec- 
ond tau decay yields one charged particle. This require- 
ment on the second tau decay greatly reduces the back- 
ground from continuum multi-hadron events. Events 
with four well-reconstructed tracks and zero net charge 
are selected, and the tracks are required to point toward 
a common region consistent with t + t~ production and 
decay. The polar angle of all four tracks in the laboratory 
frame is required to be within the calorimeter acceptance. 
Pairs of oppositely-charged tracks are ignored if their in- 
variant mass, assuming electron mass hypotheses, is less 
than 30 MeV/c 2 . Such tracks are likely to be from photon 
conversions in the traversed material. The event is di- 
vided into hemispheres in the e + e~ cm. frame using the 
plane perpendicular to the thrust axis, as calculated from 



the observed tracks and neutral energy deposits. The 
signal (3-prong) hemisphere must contain exactly three 
tracks while the other (1-prong) hemisphere must con- 
tain exactly one. Each of the charged particles found in 
the 3-prong hemisphere must be identified as a lepton or 
hadron candidate appropriate to the search channel. The 
relevant particle identification capabilities of the BABAR 
detector are described in Ref. Q. 

To further suppress backgrounds from quark pair pro- 
duction, Bhabha scattering events, and SM tau pair pro- 
duction, we apply additional selection criteria separately 
in the eight different search channels. Specific cut values 
are shown in Tab.Hl All selection criteria are optimized to 
provide the smallest expected upper limit on the branch- 
ing fraction in the background-only hypothesis. Reso- 
nant decays are selected with cuts on the invariant mass 
of the two hadrons in the 3-prong hemisphere (nihh)- 
The invariant mass of the 1-prong hemisphere (mi_ F ) 
is calculated from the charged and neutral particles in 
that hemisphere and the total missing momentum in the 
event. As the missing momentum in signal events results 
from one or more neutrinos in the 1-prong hemisphere, 
this mass is required to be near the tau mass. Back- 
ground events from quark pair production are suppressed 
with cuts on the missing transverse momentum in the 
event (p™ ss ), the scalar sum of all transverse momenta 
in the cm. frame (p^™ s ), and the number of photons in 



the 1-prong and 3-prong hemispheres (nl , 



l 3pr 



). To re- 



duce the background contribution from radiative Bhabha 
and di-muon events, the 1-prong and 3-prong momentum 
vectors must not be collinear in the cm. frame. For the 
same reason, the 1-prong track must not be identified as 
an electron for the t~ — > e~ p search. 



TABLE I: Values of the cuts on the selection variables de- 
scribed in the text. Masses are in units of GeV/c 2 , and mo- 
menta in units of GeV/c. 



Channel 


ecf> 


p(j> 


ep 


pp 


eK* pK* eK* fj,K* 


mhh min 


1.000 


1.005 


0.6 


0.6 


0.8 


0.82 


0.80 0.78 


m h h max 


1.040 


1.035 


0.92 


0.96 


1.0 


0.98 


1.04 1.00 


mi-pr min 


0.3 


0.4 


0.3 


0.3 


0.3 


0.2 


0.3 - 


mi-pr max 


2.5 


2.5 


2.5 


2.5 


2.5 


2.5 


2.5 - 


P rp lss min 


0.4 


0.3 


0.4 


0.4 


0.4 


0.4 


0.4 0.4 


p c T ms min 


0.5 








0.6 




0.3 - 


n"{ pr max 


4 


3 


3 


1 




3 


2 


rig pr max 


3 


1 


2 


1 




2 


1 



As a final discriminant, we require candidate signal 
events to have an invariant mass and total energy in 
the 3-prong hemisphere consistent with a parent tau lep- 
ton. These quantities are calculated from the measured 
track momenta, assuming lepton and hadron masses 
that correspond to the neutrinoless tau decay in each 
search channel. The energy difference is defined as 
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AE = E *cc ~ E kam. where E * C c is the total energy 
of the tracks observed in the 3-prong hemisphere and 
i?beam is the beam energy, with both quantities calcu- 
lated in the cm. frame. The mass difference is defined 
as AM = Mec — m T where Mec is calculated from a 
kinematic fit to the 3-prong track momenta with the en- 
ergy constrained to be -J ' s/2 in the cm. frame, and 
m T = 1.777 GcV/c 2 is the tau mass [12]. While the en- 
ergy constraint significantly reduces the spread of AM 
values, it also introduces a correlation between AM and 
AE, which must be taken into account when fitting dis- 
tributions in this 2-dimensional space. 

Detector resolution and radiative effects broaden the 
signal distributions in the (AM, AE) plane. Because of 
the correlation between AM and AE, the radiation of 
photons from the incoming e + e~ particles produces a tail 
at positive values of AM and negative values of AE. Ra- 
diation from the final-state leptons, which is more likely 
for electrons than muons, leads to a tail at low values 
of AE. Rectangular signal boxes (SB) in the (AM, AE) 
plane are defined separately for each search channel. As 
with previous selection criteria, the SB boundaries are 
chosen to provide the smallest expected upper limit on 
the branching fraction. The expected upper limit is esti- 
mated using only MC simulations and data events in the 
sideband region, as described below. Figure [T] shows the 
observed data in the Large Box (LB) of the (AM, AE) 
plane, along with the SB boundaries and the expected 
signal distributions. Table HT1 lists the channel-specific di- 
mensions of the SB. While a small fraction of the signal 
events lie outside the SB, the effect on the final result 
is negligible. To avoid bias, we use a blinded analysis 
procedure with the number of data events in the SB re- 
maining unknown until the selection criteria are finalized 
and all crosschecks are performed. 



TABLE II: Signal Box boundaries; AM is in units of GeV/c 2 
and AE in units of GeV. 



Mode 


ecj> 


ep 


eK* 


eK* 


n4> 


pp 


pK* 


pK* 


AMmin 


-0.02 


-0.02 


-0.02 


-0.015 


-0.008 


-0.01 


-0.01 


-0.008 


AM max 


0.015 


0.02 


0.02 


0.02 


0.01 


0.015 


0.01 


0.01 


A_B m in 


-0.13 


-0.10 


-0.15 


-0.125 


-0.09 


-0.06 


-0.08 


-0.08 


A£^ m ax 


0.10 


0.06 


0.08 


0.06 


0.06 


0.04 


0.04 


0.06 



There are three main classes of background events 
remaining after the selection criteria are applied: 
charm quark production (cc) , low-multiplicity continuum 
e + e~ — > uu/dd/ss events (uds), and SM t + t" pair 
events. The background from two-photon production 
is negligible. These three background classes have dis- 
tinctive distributions in the (AM, AE) plane. The uds 
events tend to populate the plane evenly, with a fall-off 
at positive values of AE. Events in the cc sample ex- 
hibit peaks at positive values of AM due to D and D s 
mesons, and are generally restricted to negative values 



of AE. The r + r~ background events are restricted to 
negative values of both AE and AM. 

The expected background rates in the SB are deter- 
mined by fitting a set of 2-dimensional probability den- 
sity functions (PDFs) to the observed data in the grand 
sideband (GS) region of the (AM, AE) plane. The GS 
region is defined as the LB minus the SB. The shapes 
of the PDFs are determined by fits to the (AM, AE) 
distributions of background MC samples in the LB, as 
described in Ref. [l|. The present analysis makes use of 
the same parameterization as Ref. [l[ for the AE spec- 
tra, except for the case of the cc spectrum in some search 
channels. In these cases, combinations of polynomial and 
Gaussian functions are used. The choice of PDF for the 
AM spectrum of the uds samples is the same as used in 
Ref. [lj, while the t + t~ and cc AM spectra are modeled 
with Gaussian and polynomial functions, or the Crystal 
Ball function [l6j]. All shape parameters, including a ro- 
tation angle accounting for the correlation between AE 
and AM, are determined from the fits to MC samples. 

Once the shapes of the three background PDFs are de- 
termined, an unbinned extended maximum likelihood fit 
to the data in the GS region is used to find the expected 
background count in the SB. The fits to the background 
MC samples and to data are performed separately for 
each of the eight search channels. 

We estimate the signal event selection efficiency with 
a MC simulation of lepton- flavor violating tau decays. 
Between 20% and 40% of the MC signal events pass the 
3-1 topology requirement. The efficiency for identifica- 
tion of the three final-state particles ranges from 42% for 
t~ — > fj,~K* to 82% for t~ — * e~p. The total efficiency 
for signal events to be found in the SB is shown in Ta- 
ble lllll and ranges from 4.1% to 8.0%. This efficiency in- 
cludes the branching fraction for the vector meson decay 
to charged daughters, as well as the branching fraction 
for 1-prong tau decays. 

The particle identification efficiencies and misidenti- 
fication probabilities have been measured with control 
samples both for data and MC events, as a function 
of particle momentum, polar angle, and azimuthal an- 
gle in the laboratory frame. The systematic uncertain- 
ties related to the particle identification have been esti- 
mated from the statistical uncertainty of the efficiency 
measurements and from the difference between data and 
MC efficiencies. These uncertainties range from 1.7% for 
t~ — > e~p to 9.0% for t~ — » pT p [17]. The modeling 
of the tracking efficiency and the uncertainty from the 
1-prong tau branching fraction each contribute an ad- 
ditional 1% uncertainty. Furthermore, the uncertainty 
on the intermediate branching fractions B((p,K ,K — > 
h + h~) contributes a 1% uncertainty. All other sources 
of uncertainty in the signal efficiency are found to be 
negligible, including the statistical limitations of the MC 
signal samples, modeling of radiative effects by the gen- 
erator, track momentum resolution, trigger performance, 
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FIG. 1: Observed data shown as dots in the Large Box of the (AM, AE) plane and the boundaries of the Signal Box. The 
dark and light shading indicates contours containing 50% and 90% of the selected MC signal events, respectively. 



and the choice of observables used for event selection. 

Since the background levels are extracted directly from 
the data, systematic uncertainties on the background es- 
timation are directly related to the background param- 
eterization and the fit technique used. Uncertainties re- 
lated to the fits to the background samples are estimated 
by varying the background shape parameters according 
to the covariance matrix and repeating the fits, and range 
from 3.8% to 10%. Uncertainties related to the fits for 
the background yields in the GS are estimated by vary- 
ing the yields within their errors, and range from 4.1% to 
16%. The total uncertainty on the background estimates 
is shown in Table IIIII Crosschecks of the background 
estimation are performed by comparing the number of 
events expected and observed in sideband regions imme- 
diately neighboring the SB for each search channel. No 
major discrepancies are observed. 

The number of events observed (iV b s ) and the number 
of background events expected (iVbgd) are shown in Ta- 
ble HTH The POLE calculator [lj], based on the method 
of Feldman and Cousins [l9[ , is used to place 90% CL up- 
per limits on the number of signal events (N^), which 
include uncertainties on iVbgd and on the selection effi- 
ciency (e). For the r~ — > [iT <\> search, the POLE cal- 
culation results in a two-sided interval at 90% CL for 
the number of signal events: [0.39 — 8.65]. Upper lim- 
its on the branching fractions are calculated according to 
2?ul = N^/(2eLa TT ), where the values L and a TT are 
the integrated luminosity and t + t~ cross section, respec- 
tively. The uncertainty on the product Ltr rr is 1.0%. Ta- 
blc lllll lists the upper limits on the branching fractions, as 
well as the expected upper limit B®® p , defined as the mean 
upper limit expected in the background-only hypothesis. 



TABLE III: Efficiency estimate, number of expected back- 
ground events (iVbgd), number of observed events (N ^ B ), ob- 
served upper limit at 90% CL on the number of signal events 
(Aul), ex P ec t e d branching fraction upper limit at 90% CL 
(£?cxp)i and observed branching fraction upper limit at 90% 
CL (Bu° L ). Bl% and SfjL are in units of 10" s . 



Mode 


e[%] 


iVbgd 


Nobs 


JV UL 


S 90 


R 90 


e<j> 


6.43 ±0.16 


0.68 ±0.12 





1.8 


5.0 


3.1 


pcj> 


5.18 ±0.27 


2.76 ±0.16 


6 


8.7 


8.2 


19 


ep 


7.31 ±0.18 


1.32 ±0.17 


1 


3.1 


4.9 


4.6 


fip 


4.52 ± 0.41 


2.04 ±0.19 





1.1 


8.9 


2.6 


eK* 


8.00 ±0.19 


1.65 ±0.23 


2 


4.3 


4.8 


5.9 


fxK* 


4.57 ±0.36 


1.79 ±0.21 


4 


7.1 


8.5 


17 


eK* 


7.76 ±0.18 


2.76 ±0.28 


2 


3.2 


5.4 


4.6 


pK* 


4.11 ±0.32 


1.72 ±0.17 


1 


2.7 


9.3 


7.3 



The 90% CL upper limits on the r -> I <j)/p/K*/K 
branching fractions are in the range (2.6 — 19) x 10~ 8 , 
and these limits represent improvements over the previ- 
ous experimental bounds 0] in almost all search channels. 
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